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We report the crystal structure, magnetization measurements, and band-structure calculations 
for the spin-i quantum magnet CaCu2 (SeOs^CL;. The magnetic behavior of this compound is 
well reproduced by a uniform spin-i chain model with the nearest-neighbor exchange of about 
133 K. Due to the peculiar crystal structure, spin chains run in the direction almost perpendicular 
to the structural chains. We find an exotic regime of frustrated interchain couplings owing to 
two inequivalent exchanges of 10 K each. Peculiar superexchange paths grant an opportunity to 
investigate bond-randomness effects under partial Cl-Br substitution. 
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I. INTRODUCTION 

Low-dimensional magnets remain an attractive play- 
ground to study quantum phenomenal and to under- 
stand strongly correlated electronic systems on a model 
levels Theoretical investigations disclose interesting fea- 
tures of numerous simple spin networks, such as a dia- 
mond chain^ a kagome lattice^ or a pyrochlore latticed 
The transfer of these spin lattices to real systems and 
the subsequent experimental verification of theoretical 
results are, however, rather problematic and stimulate 
extensive studies of low-dimensional (or potentially low- 
dimensional) magnetic materials. Most of these studies 
are focused on Cu +2 compounds, because the d 9 nature 
and the pronounced Jahn- Teller effect of Cu +2 ion lead 
to insulating spin-i compounds with diverse spin-lattice 
geometries. 

Aiming to find hitherto unexplored examples of low- 
dimensional spin systems, we investigate copper selenite- 
chlorides. These compounds combine two important 
structural ingredients that lead to unusual physical be- 
havior: i) Se03 selenite groups contain the lone-pair Se +4 
cation that induces exotic (and potentially polar) crystal 
structures, as in the piezoelectric ferrimagnet Cu20Se03 
(Refs. la); ii) CI atoms show strong hybridization with 
Cu 3d orbitals and mediate long-range exchange cou- 
plings, leading to highly entangled spin lattices in spin- 
tetrahedra compounds Cu2Te205X 2 (X = CI, Br) 7 - show- 
ing incommensurate magnetic order, or in the intricate 
spin-dimer system (CuCl)LaNb 2 7 (Ref. H and i) . 

Here, we present an experimental and computational 
study of CaCu2(Se03)2Cl2- X-ray diffraction, mag- 
netization measurements, and band structure calcula- 
tions are applied to elucidate crystal structure, electronic 
structure, and magnetic behavior of this compound. In 
our study, the complex structural arrangement of Cu 
polyhedra is readily disentangled by a microscopic ap- 
proach. We find the conventional x 2 — y 2 orbital ground 



state for both Cu sites, and establish a minimum mag- 
netic model of uniform spin- 1/2 chains with weak, frus- 
trated interchain couplings. 




FIG. 1. (Color online) Left: crystal structure of 
CaCu2(Se03)2Cl2- The structural chains run along [101] (not 
shown), whereas the magnetic chains range along ~ [201], as 
shown by bold (red) lines. Small unlabeled spheres denote O 
atoms. Right: local environment of Cu(l) and Cu(2). The 
magnetically active Cu(l)04 and Cu(2)02Cl2 plaquettes are 
highlighted. 



II. SYNTHESIS AND SAMPLE 
CHARACTERIZATION 

Calcium selenite CaSe03 was prepared via solution 
synthesis, as described in Ref. [PJ, The solutions of 
calcium nitrate Ca(NOs)2 (chemically pure) (6.214 g) 
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and selenous acid H 2 Se03 (98%) (4.886 g) in a mini- 
mal amount of hot distilled water were mixed. The 
ammonia 1:5 water solution was added to hx pH of 
the solution in the range 7—8. The fine white pow- 
der was obtained as a precipitate. The precipitate 
was then dried at 150 °C. According to X-ray powder 
diffraction (XRPD), the obtained powder was identi- 
fied as a hydrate CaSeOs-H^O. The hydrate was fur- 
ther calcined on a gas burner in a ceramic plate for 
30 min. The resulting product was identified as a sin- 
gle phase CaSe0 3 [space group P2 1 /n, a = 6.399(5) A, 
6 = 6.782(4)1, c = 6.682(8)A, /3 = 102.84(6)°]. 

ScC>2 was obtained from H2SCO3 by its decomposition 
under vacuum at 60 °C and the sublimation of the result- 
ing substance in a flow of anhydrous air and NO2. 

CuO (ultra pure) and CuCl 2 (Merck, >98%) 
were used. The dark-greenish powder sample of 
CaCu2(Se03)2Cl2 was obtained from a stoichiometric 
mixture of CaSeC>3, CuCk;, CuO, and Se02- The mix- 
ture (about 0.5 g total) was prepared in Ar- filled camera, 
sealed in a quartz tube, and placed into the electronically 
controlled furnace. The sample was heated from room 
temperature to 300 °C for 12 hours, exposed at 300 °C 
for 24 hours, heated up to 500 °C for 12 hours, and ex- 
posed at 500 °C for 96 hours. 

The resulting samples were single-phase, as con- 
firmed by powder x-ray diffraction (STOE STADI- 
P diffractometer, CuK a i radiation, transmission ge- 
ometry). The powder pattern was fully indexed in 
the monoclinic space group C2/c with lattice pa- 
rameters a = 12.752(3) A, 6 = 9.036(2)1, c = 6.970(1) A, 
P = 91.02(1)°. CaCu2(Se0 3 ) 2 Cl2 is rather stable in air, 
although a prolonged exposure of about 3 months led to a 
partial decomposition towards crystalline CuSe03-2H 2 
and possible amorphous products. 



III. CRYSTAL STRUCTURE 

For the structure determination, a single crystal was 
picked up from the bulk polycrystalline sample. The 
data were collected at the CAD-4 (Nonius) diffractome- 
ter (MoK a radiation) at room temperature. The analy- 
sis of systematic extinctions unambiguously pointed to 
the space group C2/c (15). The lattice parameters 
a = 12.759(3) A, 6 = 9.0450(18) A, c = 6.9770(14) A and 
/3 = 91.03(3) were refined, based on 24 well-centered re- 
flections in the angular range 12.01 o <6»<15.84°. The 
diffraction data were collected in an u — 29 mode with the 
data collection parameters listed in Table |H A semiem- 
pirical absorption correction was applied to the data 
based on ?/>-scans of seven reflections with A angles close 
to 90°. 

Positions of metal and selenium atoms were found by 
direct methods (shelxs-97).— Oxygen atoms were local- 
ized by a sequence of least-square cycles and difference 
Fourier syntheses Ap(x,y, z). The final refinement with 
anisotropic atomic displacement parameters was based 



TABLE I. Data collection and structure refinement parame- 
ters for CaCu 2 (Se0 3 )2Cl2. 
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on F 2 (shelxl-97). 11 Further information and the re- 
finement residuals are given in Table |T] Atomic coordi- 
nates and atomic displacement parameters are listed in 
Table IB^ 

The crystal structure of CaCu2(Se03)2Ci2 is shown 
in Fig. [T] It comprises two inequivalent Cu positions 
(Table [TT|- The Cu(l) atoms show a slightly distorted 
square-planar Cu(l)C>4 environment, typical for Cu +2 
oxides (Table IIII[) . In contrast, Cu(2) is six- fold coor- 
dinated, with four oxygens and two chlorines forming an 
octahedron, squeezed along Cu(2)-01 ( Table UlT]) . Al- 
though such octahedral coordination of Cu(2) is rather 
unusual, it still allows for a square-planar-like crystal- 
field splitting of 3d levels and the conventional non- 
degenerate 3c? 9 orbital ground state with the half-filled 
3d x 2_ y 2 orbital lying in the plane of a Cu(2)02Ci2 pla- 
quette (Fig.[T])ji2. This plaquette is formed by two Cu(2)- 
Cl bonds and two Cu(2)-01 bonds. The formation of the 
plaquette can be qualitatively understood in terms of dif- 
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TABLE II. Atomic coordinates and atomic displacement pa- 
rameters Ui B0 . 
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third of the trace of the orthogonalized Uij tensor. 



TABLE III. Selected interatomic distances (in A) in the 
Ca0u2(SeO3)2Cl2 structure. 
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Ca-Cl (x2) 


2.948(1) 
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Se-Ol 
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Se-03 
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ferent ionic radii for oxygen and chlorine. The larger size 
of the CI atoms makes their effect on the Cu 3d orbitals 
similar to the effect of Of with shorter distances to Cu. 
The resulting crystal-field splitting resembles that of a 
Cu04 plaquette and drives one of the atomic d orbitals 
half-filled as well as magnetic, that is a posteriori con- 
firmed by our DFT calculations (Sec. |V|) . 

In general, the formation of CUO2CI2 plaquettes is typ- 
ical for copper oxy chlorides. 9,14 However, a unique fea- 
ture of CaCu2(Se03)2Ci2 is the presence of two longer 
Cu(2)-03 bonds which look similar to the Cu(2)-Cl 
bonds in terms of interatomic distances but are essen- 
tially inactive with respect to the magnetism, as will be 
shown in Sec. IIV1 

The Cu(f)04 plaquettes and the Cu(2)04Cl2 octahe- 
dra share corners and form chains along [101]. However, 
the bridging 03 atoms do not belong to the Cu(2)02Cl2 
plaquettes, hence a simple Cu(l)-03-Cu(2) superex- 
change is unlikely. Instead, the leading exchange cou- 
plings should run via Se03 trigonal pyramids which 
join the plaquettes into a framework. The CI atoms 
shape tunnels that run along c and accommodate the 
Ca cations. Surprisingly, CaCu^SeOs^CL; bears no re- 
lation to SrCu2(Se03)2Cl2 (Ref. fl5l) owing to the smaller 
size of Ca and the high flexibility of Cu-Se-O-Cl frame- 
work. The arrangement of polyhedra does not resemble 
any known structure type either. 



IV. MAGNETIC PROPERTIES 

Magnetic susceptibility (x) of CaCu2(Se03)2Ci2 was 
measured with a Quantum Design MPMS SQUID mag- 
netometer in the temperature range 2 — 380 K in applied 
fields [i H of 0.5, 2, and 5 T. 

The x(T) dependence (Fig. [5]) shows a broad maxi- 
mum at 83 K and a pronounced increase below 30 K. The 
susceptibility maximum is a signature of quantum fluc- 
tuations (low-dimensional and/or frustrated behavior), 
while the low-temperature upturn is caused by the para- 
magnetic contribution of defects and impurities. Above 
230 K, we fit the data with the modified Curie- Weiss law 
x = Xo +C/{T+9) where xo = 6(f) x 10~ 5 emu (mol Cu)- 1 
accounts for the temperature- independent (e.g., van 
Vleck) contribution, C = 0.42(f) emu K (mol Cu) -1 is the 
Curie constant, and 9 — 93(5) K is the Weiss temperature. 
The positive 9 indicates predominant antiferromagnetic 
(AFM) interactions in the system. Using the expressions 

c= Na^s^ Meff = <?/WS(S + l) (1) 

we obtain the resulting effective magnetic moment 
/i ff = f .83(f) hb and the ^-factor g = 2.11(1), typical for 
spin-i Cu+ 2 ^ 

To fit the whole x(T) dependence, we used different 
expressions for simple low-dimensional spin models. The 
best fit was obtained with the expression for the uniform 
spin- 1 chain Xch, given by Ref. [l7| [see their Eq. (53) pa- 
rameterized with the values provided in the third column 
of Table I] . The temperature range 2 — 380 K fits to the 
validity condition of this parameterization < J <5»^ 
To account for temperature-independent and the low- 
temperature impurity contribution to x(^)> we supple- 
mented Xch with the xo term and the Curie term C lmp /T, 
respectively: 

m , gimp , N A g 2 fi% (T\ 
x(T)=Xo + ^r L + j~ Xch [j)' 

The fit yielded the intrachain exchange cou- 
pling J = f 33(f) K, the g-factor 5 = 2.11(f), 
Xo = 3(f)xfO- 5 emu(molCuT 1 f 1 £ and 
Cimp = 0.005(1) emu K (mol Cu)" 1 (about 1% of spin- \ 
impurities) . To check the applicability of the Heisenberg 
chain model to our system, we calculate x^ ax T max (7~ 2 , 
that should amount to 0.0353229(3) emuK (mol Cu)" 1 
for a Heisenberg chain system, independent of J 
(see Eq. 31 from Ref. E3). For CaCu2(Se0 3 ) 2 Cl2, 
Xoh(T max )T max g- 2 = 0.0345(8) emuK (mol Cu)" 1 devi- 
ates only by few percent from the ideal value, justifying 
the mapping onto the Heisenberg spin chain model. 

The extrinsic nature of the low-tempcrature Curie tail 
is supported by its suppression in magnetic field (Fig. [3]). 
Temperature derivative of magnetic susceptibility ex- 
hibits a kink at 6 K, which is likely a signature of antifer- 
romagnetic ordering. This issue is discussed in context 
of the microscopic spin model in Sec. I VII 
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FIG. 2. (Color online) Magnetic susceptibility (circles) and 
the fit (dashed line), see text for details. Intrinsic, i.e. Heisen- 
berg chain (solid line) and impurity (dotted line) contribu- 
tions to the fitted curve are shown. Inset: high-field magne- 
tization curve (squares) with a linear fit (dashed line). 




30 60 90 120 

Temperature (K) 

FIG. 3. (Color online) Magnetic susceptibility of 
CaCu2(SeOs)2Cl2 measured in applied fields of 0.5 T, 2T, 
and 5 T. The increase in the field leads to suppression of the 
low-temperature paramagnetic upturn. The inset shows the 
derivative of the magnetic susceptibility and a kink around 
6K, likely evidencing long-range magnetic ordering. 



We also measured the magnetization curve of 
CaCu2(Se03)2Cl2 in pulsed magnetic fields up to 60 T 
at a constant temperature of 1.5 Kji£ The linear change 
in the magnetization (see the inset of Fig. [5]) is consis- 
tent with the proposed uniform-chain behavior, since the 
accessible field range is well below the saturation field 
(jjlqHs = 188 T for J = 133 K and g = 2.11) for the param- 
eters obtained above^S 

According to Sec. IIII1 the spin chains cannot be as- 
signed to the structural chains of the corner-sharing 
Cu(l)04 plaquettes and CUO4CI2 octahedra. In Sec. \V\ 
we show that the uniform chains originate from a non- 
trivial superexchange pathway via SeC>3 groups. 



V. MICROSCOPIC MODEL 

Band structure calculations have been performed us- 
ing the full-potential code fplo9. 00-31. 21 For the ex- 
change and correlation potential, the parameterization 
of Perdew and Wang has been chosen^ For the calcu- 
lations within the local density approximation (LDA), a 
well converged /c-mesh of 10 x 10 x 12 points was used. For 
spin-polarized supercell local spin density approximation 
(LSDA)+£/ calculations, we used fc-meshes of 4x4x4 and 
4x4x2 points. Convergence of total energy with respect 
to the /c-mesh has been carefully checked. 

LDA yields a valence band (Fig.|4j top) with the band- 
width of about 8 eV, typical for cuprates^ - — The band 
is clearly split into two parts: the region between —8 and 
— 5.5eV is dominated by Se Ap and O 2p states, while the 
rest of the valence band is formed by Cu, O, and CI states. 
Non-zero density of states (DOS) at the Fermi level £f 
indicates a metallic solution in contrast to the insulating 
behavior, expected from the green sample color. This 
well-known drawback of the LDA arises from a strong un- 
derestimation of correlations which are intrinsic for the 
3d 9 electronic configuration (Cu 2+ ) and drive the system 
into the insulating regime.— 

The states relevant for magnetism are confined to the 
vicinity of ef- In most cuprates, these are the antibond- 
ing Cu 3d x 2_ y 2 and O 2p a states (in the local coordinate 
system) , typically well-separated from the rest of the va- 
lence band£2&2!28 However, CaCu 2 (Se0 3 )2Cl2 lacks a 
separated band complex around ep. This reflects the 
octahedral coordination for Cu(2) and makes a detailed 
analysis of the magnetically active orbitals necessary. 

To evaluate the relevant states, we project the DOS 
onto a set of local orbitals. This way, we find the 
dominant 3d x 2_ y 2 contribution to the Cu(l) DOS at E-p 
(Fig.lU left bottom), as for almost all undoped cuprates. 
For the Cu(2) atom, the situation is less trivial, since the 
local environment of this atom implies two short Cu(2)-0 
bonds as well as four long (two Cu(2)-0 and two Cu(2)- 
Cl) almost equidistant bonds, making the choice of the 
local coordinate system ambiguous. However, the analy- 
sis of local DOS for different situations readily yields the 
correct choice of the local axes and evidences that the 
two short Cu(2)-0 and two Cu(2)-Cl bonds form a pla- 
quette with the Cu 3d x 2_ y 2 magnetically active orbital 
(Fig. [5]). The local DOS of this orbital clearly dominates 
the states at £f (Fig- IH right bottom) and confirms our 
empirical considerations presented in Sec. IIIII 

Since the magnetism of CaCu2(Se03)2Cl2 is confined 
to Cu 3d x 2_ y 2 orbitals, they can be used as a minimal 
basis for an effective tight-binding (TB) model. The to- 
tal number of states (four) in the model corresponds to 
the number of magnetic Cu atoms in a unit cell: two 
Cu(l) and two Cu(2). To parameterize the model, we 
use the Wannier functions (WFs) technique, which yields 
numerical values for the leading hoppings tj (transfer in- 
tegrals). This way, we obtain a perfect fit to the LDA 
bands (Fig. 2J bottom). 
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FIG. 4. (Color online) Top: total (shaded) and atomic- 
resolved (lines) density of states (DOS) for CaCu2(Se03)2Cl2. 
Bottom center: LDA band structure (circles) and the WF fit 
(dashed line). Bottom left and right: orbital-resolved density 
of states for Cu(l) and Cu(2). Solid [dashed] lines denote total 
DOS for Cu(l) [Cu(2)]. Shaded regions show the 3d x 2_ y 2 con- 



tribution. O 2p c 
(dotted) lines. 



(CI 3p CT ) DOS are shown with dash-dotted 



TABLE IV. Interatomic Cu-Cu distances (d, in A), transfer 
integrals (i;, in meV), as well as antiferromagnetic (AFM) and 
ferromagnetic (FM) contributions to the total exchange inte- 
grals (Ji, in K) for the leading couplings in CaCu2(Se03)2Ci2. 
For notation of the paths see Fig.[7]( J nn stands for the nearest- 
neighbor coupling along the structural chains). 
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FIG. 5. (Color online) Band weights for the magnetically 
active Cu 3d x 2_ y 2 and other Cu 3d orbitals. 



Within the effective one-orbital TB model, we find 
three relevant couplings (Table HV*)) : t\ (Figs. [5] and [7]) 
running along Cu(l)-Cu(2) chains almost parallel to the 
[201] direction, the short-range interchain coupling t\ c \ 
which connects Cu(l) atoms, as well as the long-range 
interchain coupling i; c2 connecting Cu(2) atoms along 
[101]. The clearly dominant t\ amounts to 139 meV, 
while iid and t- lC 2 are found to be 47 and 30meV, re- 
spectively. Due to the particular orientation of magnet- 
ically active orbitals, the hopping t nn along the "struc- 
tural chains" is apparently small, does not affect the mag- 
netic ground state, as will be shown later. 

To restore the insulating ground state, we map our 
TB model onto a Hubbard model with an effective on- 



FIG. 6. (Color online) Wannier functions marking the Ji 
superexchange path. Cu(l)04 (left) and Cu(2)02Ci2 (right) 
plaquettes are filled. The Se03 pyramid is visualized by Se-O 
bonds (lines). The projection is the same as in Fig. [Tj (left). 



site Coulomb repulsion U c s. Since the strongly cor- 
related limit (U c tf^>ti) and the half-filling are well- 
justified for undoped cuprates, the low-energy (mag- 
netic) excitations of the Hubbard model can be well 
described by a Heisenberg model. Within the second- 
order perturbation theory, AFM exchange couplings 
are expressed as J AFM = Atf /U c s. 29 Adopting a typi- 
cal value U eS = 4.5eV£&2± we obtain J AFM = 200K, 
jA™ = 25 K and J AFM = 10 K. 

The one-orbital approach yields only the AFM part 
jAfm £ total exchange Ji. However, particular ge- 
ometrical configurations (Cu-O-Cu angle close to 90°, 
edge-sharing of CUO4 plaquettes etc.) often lead to siz- 
able ferromagnetic (FM) contributions to the total ex- 
change, which can even outgrow </ AFM , resulting in a FM 
(negative) Ji. Based on structural considerations only, an 
appreciable FM contribution might be expected for the 
short-range coupling J lc \ and the nearest-neighbor cou- 
pling J nn , whereas J\ and J 1C 2 are rather long-range, and 
their FM contributions should be small. To challenge 
this conjecture, we perform supcrccll LSDA+L 7 calcula- 
tions. This method gives access to the total exchange Ji, 
being the sum of the AFM and FM contributions, J AFM 
and Jf M , respectively. Combining the LSDA+C/ results 
with J AFM estimates from the model approach described 
above, the FM contributions jf M can be evaluated. 



6 



Recent thorough studies on low-dimensional cuprates 
give evidence that quantitative magnetic models based 
on the results of LSDA+17 calculations are, in addi- 
tion to the dependence on the Coulomb repulsion Ud, 
also dependent on the double-counting correction (DCC) 
scheme . 25 ' 26 The influence of these two parameters is par- 
ticularly large for systems with sizable Jf M - Yet, for 
cuprates with structurally isolated plaquettes, as in the 
case of CaCu2(Se03)2Ci2, the around-mean-field (AMF) 
DCC with the Ud value 6.5 ± 0.5 eV typically yields accu- 
rate (±10%) estimates for the leading couplings. Adopt- 
ing the AMF DCC and U d = 6.5 eV, we obtain J x = 145 K, 
in excellent agreement with the experimental J ~ 133 K 
from the fit to the magnetic susceptibility. In accordance 
with our expectations, the long-range interchain cou- 
pling Jj C 2 has a tiny FM contribution only, while for the 
short-range coupling J- lc i the FM contribution reaches 
Jf c f = -15K (see Table E}. At first glance, the rather 
large jf M ~ — 55 K may look surprising. However, in the 
related quasi-lD system CuSe20s, a similar intrachain 
coupling runs via two corner-sharing Se03 pyramids and 
the FM contribution jf M to this coupling amounts to 
even — lOOK^i Therefore, such high FM contributions 
are likely intrinsic to the superexchange realized via Se03 
pyramids. Understanding the underlying mechanism of 
this complex superexchange deserves a separate study 
and lies beyond the scope of the present investigation. 
Presently, we note that the superexchange results from 
the overlap of oxygen orbitals, while Se states have only 
a minor contribution to the WFs. 

The last remark concerns the short-range coupling be- 
tween the nearest-neighbor Cu(l) and Cu(2) atoms that 
form the structural chains. The WFs analysis yielded 
a negligible t associated with this coupling path. Still, 
the respective interatomic distance (3.84 A) is relatively 
small, which could give rise to a FM coupling. Therefore, 
we evaluated the total exchange J nn using the LSDA+Z7 
calculations. The resulting exchange of 4K is in excel- 
lent agreement with the TB estimate (4K), evidencing a 
negligible FM contribution and justifying our restriction 
to the three couplings Ji, J; c i, and J lC 2 for a minimum 
model. 



VI. DISCUSSION AND SUMMARY 

The spin model of CaCu2(SeC>3)2Cl2 is depicted in 
Fig. [7] Its main element are Cu(l)-Cu(2) chains running 
almost parallel to the [201] direction, which is different 
from the structural chains. The spin-^ chains are cou- 
pled by two inequivalent exchange interactions: J; c i is 
short-range and links the Cu(l) atoms of the two neigh- 
boring chains, while the long-range J lC 2 bridges the Cu(2) 
atoms of the fourth- neighbor chains. Another difference 
between J lc \ and J lC 2 is that the former is responsible 
for a 3D coupling [connects Cu(l) atoms belonging to 
different layers, see Fig. [7J, whereas the latter is con- 
fined to the ac plane. Either of inter-chain couplings 




FIG. 7. (Color online) Spin model for CaCu 2 (Se0 3 )2Cl2. 
Filled and empty circles show the Cu(l) and Cu(2) positions. 
Bold lines and circles denote couplings (see legend for no- 
tation; the same notation applies for Table IIV[) in the front 
plane, whereas gray circles and shaded lines correspond to the 
atoms lying in the back plane. The planes are connected by 
Jid couplings only. A closed loop (bold line) having an odd 
number of AFM couplings indicates that the spin model is 
frustrated. The unit cell is depicted by the gray rectangle. 



alone, J- lc i or J; C 2, leads to a 3D or 2D non- frustrated 
model, respectively. However, the combination of the 
two inter-chain couplings gives rise to magnetic frustra- 
tion, evidenced by an odd number of AFM bonds along 
the closed, hourglass-shaped loop shown in Fig. [7J 

In general, most of the quasi-one-dimensional cuprates 
order AFM at low temperatures, while strong quantum 
fluctuations drastically reduce the value of the ordered 
moment compared to the classical value of 1/LtB-— The 
AFM transition temperature can vary in a wide range, 
since it depends not only on the leading intrachain cou- 
pling Ji , but also on the absolute values and the topology 
of the interchain couplings. To emphasize the huge im- 
pact of the interchain coupling regime onto the transition 
temperature, we mention here two quasi-one-dimensional 
systems with similar intrachain coupling J\ of 150- 
200 K, but different interchain coupling regimes: a frus- 
trated interchain coupling in Sr2Cu(P04)2 leads to a very 
low ordering temperature of 85 mK (T/v/J = 4.5 xl0~ 4 ) 
fRef. l3ll) despite an interchain coupling of 3 KM On the 
contrary, a sizable and non-frustrated interchain coupling 
Ji c ~20K in CuSe205 results in a long-range AFM or- 
dering at rather high T N = 17 K (T N /J= 0.1)^ 

Since CaCu2(Se03)2Cl2 exhibits a similar energy scale 
(Ji ~ 133 K), it is natural to compare this compound to 
the aforementioned systems. The major difference here 
is the presence of two types of interchain couplings, J- lc i 
and Jic2, which form a 3D spin model, in contrast to 
Sr2Cu(P04)2 and CuSe2 05, where the leading interchain 
couplings are confined to 2D, while the coupling along 
the third direction J± is substantially smaller. This ar- 
gument favors higher Tjv in CaCu2(Se03)2Cl2. On the 
other hand, the interchain couplings in CaCu2(Sc03)2Cl2 
are frustrated, which certainly inhibits the magnetic or- 
dering and can considerably lower the Tjy. We expect the 
combination of the 3D coupling regime and the frustra- 
tion to result in a moderate Tjy of CaCu2(Se03)2Cl2, 
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comparable to that of CuSe20s. Indeed, the kink of 
magnetic susceptibility at 6 K in Fig. [3] fits well to the 
energy scale of the anticipated long-range ordering in 
CaCu 2 (Se0 3 )2Cl2 (compare to Tjv = 17K in CuSe 2 5 
with non- frustrated interchain couplings). 

In general, magnetic specific heat C™ agn data provide 
an independent estimate for the leading exchange cou- 
pling and are sensitive to the long-range magnetic or- 
dering. However, the measured specific heat contains, 
apart from C™ agn , also a phonon contribution. The en- 
ergy scale of magnetic interactions in CaCu2(Se03)2Cl2 
gives rise to a maximum in C™ asn (T) at rather high 
T« (0.48 J) = 64K (Eq. 39 in Ref.[l|). At this temper- 
ature, the phonon specific heat strongly dominates over 
C™ agn , impeding an accurate disentanglement of mag- 
netic and phonon contributions. Moreover, the large 
value of J leads to only a small amount of magnetic 
entropy, which could be released at the transition tem- 
perature. Thus, the resulting magnetic specific heat 
anomaly would be rather small^- Therefore, consider- 
ing large quantum fluctuations that substantially lower 
the ordered magnetic moment,— the method of choice 
are muon spin resonance (/iSR) experiments that should 
be carried out in future to verify magnetic ordering in 
CaCu2(SeC>3)2Cl2. For instance, long-range magnetic or- 
dering in the square-lattice compounds Cu(Pz)2 (0104)2 
and [Cu(Pz) 2 (HF 2 )]BF4 was only revealed by /iSR, while 
the heat capacity data lack any signatures of transition 
anomalies^ 

Magnetic frustration is one of the leading mechanisms 
that give rise to complex magnetic structures. It is there- 
fore interesting to address the nature of the anticipated 
magnetically ordered ground state of CaCu2(Se03)2Cl2. 
First, we consider the low-energy sector of a classical 
Heisenberg model on finite lattices of 16 coupled chains. 
For each chain, we impose a condition of the ideal an- 
tiferromagnetic arrangement of the neighboring spins. 34 
To keep the problem computationally feasible, we first 
restrict ourselves to collinear spin arrangements. The 
magnetic ground state is evalauted as a state with min- 
imal energy. Adopting the ratios of the leading ex- 
change couplings from our LSDA+[7 calculations (Ta- 
ble IIV1 last column) , we arrive at an AFM ground state, 
with the magnetic unit cell doubled along a and quadru- 
pled along c with respect to the crystallographic unit cell, 
i.e. the propagation vector is (h,0, 4). To understand 
the particular way the frustration is lifted, we analyse 
which couplings are satisfied, by considering the prod- 
ucts — A[Si • Sj](Jij/\Jij\) for all spin pairs in a unit 
cell. For collinear configurations, such product amounts 
either to 1 (a satisfied coupling) or to —1 (an unsatisfied 
coupling). For a certain type of exchange coupling, the 
sum of such products can be divided by the total number 
of couplings of this type in the unit cell (note the multi- 
plicities different from 1). This way, we can estimate the 
fraction of satisfied couplings. Such analysis yields that 
100% of Ji and J ic2 , but only 75% of J icl couplings are 
satisfied in the proposed (5, 0, |) ground state. 



Taking into account the restriction to collinear states, 
it is worth to address the stability of this ground state 
using alternative techniques. Thus, we use a classical 
Monte-Carlo code^ from the ALPS simulations package^ 
and calculate diagonal spin correlations (SfSJ), where 
i and j are spins coupled by a particular magnetic ex- 
change. This way, we obtain -0.24452(1), -0.16716(1) 
and —0.23056(3) for the Ji, J; c i and J ic2 couplings, 
respectively. These numbers should be compared to 
{S*Sj) — —0.25 for a perfect antiferromagnetic arrange- 
ment. Despite the small deviations from this ideal num- 
ber, the spins coupled by J\ and J- 1C 2 can be regarded as 
antiferromagnetically arranged, corroborating our clas- 
sical energy minimization result. On the contrary, the 
value for J lc \ is substantially smaller, yielding the average 
angle of arccos (0.16716/S' 2 ) « 48 between the respec- 
tive spins. The resulting angle is very close to 7r/4, hence 
spins in the fourth-neighbor chains are almost antiparal- 
lcl to each other (the angle amounts to 71"). This is in 
accord with the almost antiparallel arrangement of spins 
coupled by J; C 2 (coupling between the fourth-neighbor 
chains) . 

In the classical model, the exotic regime of frustrated 
interchain couplings leads to a rather complex magnetic 
ordering in CaCu2(Se03)2Ci2: the classical energy min- 
imization yields the collinear (|,0, 4) state, while the 
classical Monte Carlo simulations are in favor of a non- 
collinear magnetic ground state. These two ground states 
differ only by the mutual arrangement of spins coupled 
by Jici- 

Since for quasi-one-dimensional systems quantum fluc- 
tuations are crucial, the respective quantum model 
should be addressed. However, the study of a mag- 
netic ordering for a three-dimensional frustrated quan- 
tum magnet is a challenging task, since standard meth- 
ods, such as exact diagonalization, quantum Monte- 
Carlo, and the density matrix renormalization group 
technique, are either not applicable or do not account for 
the thermodynamic limit. Moreover, CaCu2(Se03)2Cl2 
features a non-negligible magnetic impurity contribution, 
as evidenced by the low-temperature upturn in x{T) 
(Fig. [2]). At low temperatures, these impurities can 
give rise to strong internal fields, 37 and possibly affect 
the ground state. Therefore, the magnetic ordering in 
CaCu 2 (Se03)2Cl2 deserves additional investigation us- 
ing alternative techniques, both from the experimental 
as well as the theoretical side. 

In contrast to well-known uniform-chain systems, 
CaCu2(Se03)2Cl2 shows a complex crystal structure 
with two Cu positions revealing an apparently differ- 
ent local environment. Our DFT calculations evalu- 
ate the 3d x 2_ y 2 magnetically active orbitals lying within 
the Cu(l)0 4 and Cu(2)0 2 Cl 2 plaquettes. Although the 
Cu(2)-03 bond lengths are similar to those of Cu(2)-Cl 
bonds, the symmetry of the magnetic orbitals renders the 
nearest-neighbor superexchange J nn path (along Cu(2)- 
03 bonds) essentially inactive (Table I IV I) . The spin 
chains run in a different direction which is dictated by 
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the orbital state of Cu and bythe suitable overlap of the 
ligand orbitals in the SeC>3 groups. This nontrivial situ- 
ation is very typical for spin- 1/2 systems where spin lat- 
tices are essentially decoupled from the low-dimensional 
features of the crystal structure: recall, for example, 
(CuCl)LaNb 2 7 £ BiCu 2 P0 6 ^ Cu 2 (P0 3 ) 2 CH 2 ,.22 and 
CuTe 2 5 42 

Although the intricate regime of the interchain cou- 
plings in CaCu 2 (Se03) 2 Cl 2 complicates theoretical stud- 
ies, the nontrivial structural organization of this com- 
pound also has an important advantage. The CI and Br 
atoms are known to be easily substitutable owing to their 
similar chemical nature. The substitution is commonly 
used to create bond randomness and to access the ex- 
otic behavior of partially disordered spin systems 4^ If 
the CI atoms take part in the superexchange, the chem- 
ical substitution inevitably changes the geometry of the 
superexchange pathways, and immediately leads to dra- 
matic changes in the spin system. In CaCu 2 (Se03) 2 Cl 2 , 
the CI atoms lie away from the leading superexchange 
pathway (Fig. [6]) , and a moderate alteration of the mag- 
netism should be expected. Partial CI/Br substitution 
will basically modify the relevant microscopic parame- 
ters (such as the crystal-field splitting) without changing 
the superexchange geometry. 

In summary, we have investigated the crystal struc- 
ture, electronic structure, and magnetic behavior of 
CaCu 2 (Se03) 2 Cl 2 . The compound comprises two Cu 
sites with essentially different local environment, but the 



same magnetically active orbital of local x 2 — y 2 sym- 
metry. A peculiar arrangement of magnetic plaquettes 
makes CaCu 2 (Se03) 2 Cl 2 a good realization of the spin-i 
antiferromagnetic Heisenberg chain model with an intra- 
chain exchange coupling of ~133 K and frustrated in- 
terchain couplings realized via two inequivalent superex- 
change paths. A kink in the magnetic susceptibility at 
6 K hints at long-range magnetic ordering, which is sub- 
ject to future experimental verification. The good poten- 
tial for a partial substitution of CI by Br atoms allows 
to look at the material from a different point of view. In 
particular, CI atoms located close to but not directly on 
the leading superexchange path make CaCu 2 (Se03) 2 Cl 2 
a promising model system to study the bond-randomness 
effects — like glass formation — in low-dimensional mag- 
nets. 
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